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Abstract 

A new combined NLO QCD analysis of the polarized inclusive and 
semi-inclusive DIS data is presented. In contrast to previous combined 
analyses, the 1/Q 2 terms (kinematic - target mass corrections, and 
dynamic - higher twist corrections) in the expression for the nucleon 
spin structure function g\ are taken into account. The new COMPASS 
data are included in the analysis. The impact of the semi-inclusive 
(SIDIS) data on the polarized parton densities (PDFs) and on the 
higher twist corrections is demonstrated. The controversial behavior 
of the polarized strange quark density obtained from the fit to the 
DIS data alone, and a combined analysis of DIS and SIDIS data is 
discussed. 



PACS numbers: 13.60.Hb, 12.38.-t, 14.20.Dh 



1 Introduction 



Experiments on polarized inclusive deep inelastic lepton-hadron scattering 
(DIS), reactions of the type I + p — > I' + X with both polarized lepton and 
hadron, because of the non-existence of neutrino data, can only, in principle 
yield information on the sum of quark and antiquark parton densities i.e. 
information on the polarized densities Au + Au, Ad + Ad, As + As and AG. 

Information about the antiquark densities Au, Ad and the separate As 
and As strange densities thus has to be extracted from other reactions, 
notably polarized semi-inclusive lepton-hadron reactions (SIDIS) I + p — > 
I' + h + X, where h is a detected hadron in the final state, or from semi- 
inclusive hadronic scattering (SIHS) like p + p — > h + X, involving polarized 
protons, and only possible at the RHIC collider at Brookhaven National 
Laboratory. 

In contrast to the situation in unpolarized DIS, a large portion of the 
most accurate data on polarized DIS lie in a kinematical region where Target 
Mass Corrections (TMC) of order M 2 /Q 2 (whose form is exactly known), and 
dynamical Higher Twist (HT) corrections of order Aq CD /Q 2 are important 
[H|2]. We have thus included such terms in our description of the DIS data. 
However, for the SIDIS data, we do not know the analogous results at present, 
so do not include such terms. As it happens almost all the SIDIS data 
we utilize are in kinematic regions where such corrections should not be 
important. 

In this talk we present the results of our combined NLO QCD analysis 
of polarized inclusive and semi-inclusive DIS data [3]. The new COMPASS 
data [U-Ej are also included in the fit. In the calculations of the semi- inclusive 
asymmetries A\ N {x, z,Q 2 ), the NLO MRST'02 PDFs [7j and the NLO DSS 
fragmentation functions (FFs) [8] were used for the unpolarized parton den- 
sities and the fragmentation functions, respectively. The new results for the 
polarized PDFs are compared to both the LSS'06 PDFs [TJ, obtained from 
the fit to the inclusive DIS data alone, and to those obtained from the first 
global analysis performed by DSSV group [9]. 

2 Results of analysis 

In this Section we present the numerical results of our global NLO QCD 
fit to the world inclusive and semi-inclusive DIS data (for references to the 
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data sets see our paper [3]). The data used (841 experimental points for DIS 
and 202 experimental points for SIDIS) cover the following kinematic regions: 
{0.005 < x < 0.75, 1 < Q 2 < 62 GeV 2 } for DIS and {0.005 < x < 0.48, 1 < 
Q 2 < 60 GeV 2 } for SIDIS processes. The statistical and systematic errors are 
added in quadrature and the uncertainties of the polarized PDFs presented 
correspond to Ax 2 = 1. A good description of the data is achieved for both 
the inclusive (XjvrP = 0-85) and semi-inclusive (xTv r p = 0-90) processes (NrP is 
the number of corresponding experimental points). The total value of Xdf 
is 0.88. The quality of the fit to the data is demonstrated in Fig. 1 in [3]. 

2.1 The role of semi-inclusive DIS data in determin- 
ing the polarized sea quark densities: Controversy 
about strange quark polarization 

Due to SIDIS data a flavor decomposition of the polarized sea is achieved 
and the light anti-quark polarized densities Au(x) and Ad(x) are determined 
without any additional assumptions. While Ad(x) is negative for any x in 
the measured x region, Au(x) is a positive, passes zero around x = 0.2 and 
becomes negative for large x. Sign-changing solutions are also found for the 
polarized strange sea As(x) and gluon AG(x) densities. The sign-changing 
behavior for AG(x) is not surprising since it was already found from the 
analysis of the inclusive DIS data alone On the other hand, on the 
basis of results from all published analyses of inclusive DIS, we consider the 
sign-changing solution for As(x) quite puzzling. The central values of the 
sea quark and gluon polarized densities together with their error bands are 
presented and compared to those of DSSV (dashed curves) in Fig. 1. 

Our LSS'06 PDFs (dot curves) [I] obtained from the NLO QCD analysis 
of the world inclusive DIS data are also presented in Fig. 1. While the 
light anti-quark polarized densities Au(x) and Ad(x) cannot be, in principle, 
determined from polarized inclusive DIS data, the sum (As + As)(x,Q 2 ) is 
well determined and all the NLO QCD analyses yield for this sum a negative 
value for any x in the measured region (for example, see Refs. [HEED])- 111 
these analyses, however, a term like (1 + 72;), which would permit a sign- 
change, was not included in the input parametrization of (As + As)(x, Ql)/2. 
We therefore re-analysed the world polarized inclusive DIS data using such 
a term in the input strange sea quark density 

x(As + As)(x, Ql)/2 = Ax a (l - xf{\ + 7 x). (1) 
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Figure 1. Our NLO sea quarks and gluon polarized PDFs at Q 2 = 2.5 GeV 2 in 
the MS scheme. For comparison the DSSV PDFs [9] are also presented. 



Our preliminary results confirm the previous ones, namely, that (As + 
As)(x, Q 2 )/2 is negative in the measured (x, Q 2 ) region. So, the behaviour 
of the polarized strange quark density remains controversial. Note that in 
the presence of SIDIS data As and As can, in principle, be separately deter- 
mined, as was done recently by the COMPASS Collaboration, where it was 
shown [TT] that there is no significant difference between As(x) and As(x) 
in the x-range covered by their inclusive and semi-inclusive DIS data. How- 
ever, the errors of the extracted values of the difference x(As(x) — As(x)) are 
rather large to allow us to conclude that the assumption As(x) = As(x) used 
in our's and the DSSV analyses is correct. So, if it is not correct, it might 
possibly be the cause that (As + As)(x,Q 2 )/2 densities obtained from the 
analyses of inclusive DIS data and combined inclusive and semi-inclusive DIS 
data sets, respectively, are in contradiction. However, at first glance, it looks 
as if the difference between As and As would have to be quite significant and 
might contradict the COMPASS results. Perhaps a more important issue is 
the sensitivity of the results to the form of the fragmentation functions. An 
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analysis by the COMPASS group [3] demonstrated that the determination of 
As(x) strongly depends on the set of the fragmentation functions used in the 
analysis and that the DSS FFs are crucially responsible for the unexpected 
behavior of As(x) obtained from the combined analysis. 



0.03 r 

0.02 - 

0.01 ■ 

0.00 , 

-0.01 ■ 

-0.02 ■ 

-0.03 - 



- LSS06(DIS) 
DIS+SIDIS(only A* J 
DIS+SIDIS(without A"J 

-LSS10(DIS+SIDIS) 



Q = 1 GeV 



x(as+As)/2 



xAG 



Q = 2.5 GeV 




Figure 2. Dependence of the form of the polarized strange quark density on 
different fits (left). Comparison between the positive and sign-changing gluon 
densities (right). 



In order to understand better the issue of the sign-changing behavior of 
As(x) in the case of the combined fit to the inclusive DIS and SIDIS data a 
more detailed analysis has been done. In this analysis the very recent COM- 
PASS data on the asymmetries A^ + p { ~\ A 1 / ) { ~ ) for charged pions and kaons 
produced on a proton target |TT] have been included. First, we have per- 
formed a fit to the DIS and SIDIS data including only the data on the pion 
A\ N asymmetries. Note that in this case only the sum x(As + As)(x, Qq) can 
be determined from the data because of the reasonable assumption = 
used for all the sets of the fragmentation functions. Second, we fitted the in- 
clusive DIS and SIDIS data excluding the data on the pion A\ N asymmetries. 
The results on x(As + As)/ 2 are illustrated in Fig. 2 (left) and compared 
to those obtained from DIS (dot curve) and the combined DIS and SIDIS 
(solid curve) analyses. Note that for the fragmentation functions the DSS 
set was used. As seen from Fig. 2 (left), in the presence only of the AJ N 
data x(As + As)/ 2 (dashed curve) is still negative in the measured x re- 
gion. The exclusion of these data from the full set of SIDIS data leads to 
a sign-changing behavior of x(As + As)/ 2 (dash dot curve). Note that in 
the later case the assumption As(x) = As(x) is used because the accuracy 
of the present SIDIS data is not enough to separate the strange quark and 
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anti-quark polarized densities. One can conclude from this study that the 
main reason As(x) to change a sign are the kaon data and the kaon FFs 
which are less known and very different for the different sets of fragmenta- 
tion functions. So, the study of the sensitivity of As(x) to the different kaon 
FFs used in the analysis is one of the key points we plan to investigate in 
the future. 

In Fig. 3 we present our results for the polarized Au(x) and Ad(x) 
densities at Q 2 = 2.5 GeV 2 , which are consistent with those obtained by 
DSSV (dashed blue curves). As expected, the SIDIS data do not influence 
essentially the sums (Au(x) + Au(x)) and (Ad(x) + Ad(x)) already well 
determined from the analysis of the inclusive DIS data. This fact is illustrated 
in Fig. 2 where our results from the combined analysis are compared with 
our LSS'06 PDFs. 
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Figure 3. Our NLO Am, Ad, (Au + An) and (Ad+Ad) polarized parton densities 
at Q 2 = 2.5 GeV 2 . DSSV [9] as well as LSS'06 [I] results for the corresponding 
densities are presented too. 
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2.2 The sign of gluon polarization 



We have found that the combined NLO QCD analysis of the present polarized 
inclusive DIS and SIDIS data cannot rule out the solution with a positive 
gluon polarization. The values of \ 2 /DF corresponding to the fits with 
sign-changing xAG(x, Q 2 ) and positive xAG(x, Q 2 ) are practically the same: 
X 2 /DF(node xAG) = 0.883 and x 2 /DF(xAG > 0) = 0.888, and the data 
cannot distinguish between these two solutions (see Fig. 2 (right)). The sea 
quark densities obtained in the fits with positive and sign-changing xAG(x) 
are almost identical. Note that the extracted HT values corresponding to 
both fits are also effectively identical. As a result, one can conclude that 
including the SIDIS data in the QCD analysis does not help to constrain 
better the polarized gluon density. 




Figure 4. Comparison between the experimental data and NLO (MS) curves for 
the ratio AG{x)/G{x) at Q 2 = 3 GeV 2 (left - high p t pairs) and Q 2 = 13 GeV 2 
(right - open charm) corresponding to positive and sign-changing xAG. Error 
bars represent the total (statistical and systematic) errors. The horizontal bar on 
each point shows the x-range of the measurement. The NLO DSSV [9] and AAC 
(the last Ref. in [10]) curves on AG(x)/G(x) are also presented. 



In Fig. 4 the ratio AG(x)/G(x) calculated for both the sign-changing 
and positive solutions for AG(x) obtained in our NLO QCD analysis is com- 
pared with the directly measured values of AG/G obtained from a quasi-real 
photoproduction of high p t hadron pairs [T2TTI4"] . and from the open charm 
production [15] measurements. For the unpolarized gluon density G(x) in 
the ratio above we have used that of the NLO MRST'02 [7]. The theoretical 
curves are given for fi 2 = 3 GeV 2 (high p t hadron pairs) and \j? = 13 GeV 2 
(open charm). As seen from Fig. 4, both solutions for the polarized gluon 
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density are well consistent with the experimental values of AG/G. It should 
be noted, however, that in the extraction of AG/G by the experiments a LO 
QCD treatment has been used. A NLO extraction of the measured values is 
needed in order for this comparison to be quite correct. In conclusion, the 
magnitude of the gluon density xAG(x) obtained from our combined NLO 
QCD analysis of inclusive and semi-inclusive DIS data and independently, 
from the photon-gluon fusion processes, is small in the region x ~ 0.08 — 0.2. 

When our combined analysis was finished, the COMPASS Collaboration 
reported the first data on the asymmetries Al + p { ~ ) , A 1 ^~' > for charged pions 
and kaons produced on a proton target [II] . As seen in Fig. 5, our predictions 
for these asymmetries are in a very good agreement with the data at measured 
x and Q 2 . 




Figure 5. Our predictions for the COMPASS asymmetries for charged pions and 
kaons produced on a proton target. 
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2.3 High twist effects 



In contrast to other combined analyses of the inclusive and semi-inclusive 
DIS data, we take into account the target mass and higher twist corrections 
in a the DIS sector. The values of the HT corrections to g\ extracted from 
the data in this analysis are shown in Fig. 6. Compared to the HT(LSS'06) 
corrections obtained in our analysis of the inclusive DIS data alone [I] the 
values of the HT corrections for the proton target are practically not changed, 
while the central values of HT corrections for the neutron target are smaller 
in the region x < 0.2, but in agreement with HT^ n )(LSS'06) within the errors, 
excepting the x region around x — 0.1 We consider the tendency of the HT^ 
corrections to be smaller in the region x < 0.2 to be a result of the new 
behavior of As(x) i.e. positive for x > 0.03. The positive contribution in g\ 
from As(x) should be compensated by a less positive HT^ contribution in 
this region. Since the biggest difference between the values of As(x)(dis+sidis) 
and As(x)dis is m the region x ~ 0.1 (see Fig. 1) this effect is biggest in this 
x region. The impact of As(x) on HT corrections is visible mainly for the 
neutron target because the contribution of As(x) in g\ is relatively larger 
than that in g\. 
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Figure 6. Impact of SIDIS data on HT effects for proton and neutron targets. 



Note that our results on the HT corrections to the nucleon spin struc- 
ture function gi(x,Q 2 ) are in a good agreement with the phenomenological 
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study of their first moments [ISJdT], as wen as with the QCD sum rule esti- 
mates [IB], the large N c limit in QCD [TH] and the instanton model [T9"|l2"0"] 
predictions. 

2.4 The spin sum rule 

Let us finally discuss the present status of the proton spin sum rule. Using 
the values for A£(Q 2 ) and AG(Q 2 ) at Q 2 = 4 GeV 2 , the first moments of 
the quark singlet AE(x, Q 2 ) and gluon AG(x, Q 2 ) densities, obtained in our 
analysis one finds for the spin of the proton: 



In Eq. © L Z (Q 2 ) is the sum of the angular orbital momenta of the quarks 
and gluons. Although the central values of the quark-gluon contribution in 
(j2J) are very different in the two cases, in view of the large uncertainty coming 
mainly from the gluons, one cannot yet come to a definite conclusion about 
the contribution of the orbital angular momentum to the total spin of the 
proton. 

3 Summary 

A new combined NLO QCD analysis of the polarized inclusive and semi- 
inclusive DIS data is presented. In contrast to previous combined analyses, 
the 1/Q 2 terms (kinematic - target mass corrections, and dynamic - higher 
twist corrections) to the nucleon spin structure function g 1 are taken into ac- 
count. The new results for the PDFs are compared to both the LSS'06 PDFs 
obtained from a fit to the inclusive DIS data alone, and to those obtained 
from the DSSV global analysis. The role of the semi-inclusive data in deter- 
mining the polarized sea quarks is discussed. Due to SIDIS data Au(x, Q 2 ) 
and Ad(x, Q 2 ), as well Au(x,Q 2 ) and Ad(x,Q 2 ) are determined without 
additional assumptions about the light sea quarks. The SIDIS data, anal- 
ysed under the assumption As(x, Q 2 ) = As(x,Q 2 ), imposes a sign-changing 
As(x,Q 2 ), as in the DSSV analysis, but our values are smaller in magni- 
tude, less negative at x < 0.03 and less positive for x > 0.03. Note that 




±AZ(Q 2 ) + AG(Q 2 ) + L Z (Q 2 ) 

-0.21 ± 0.46 + L Z (Q 2 ) (node AG), 
0.42 ± 0.19 + L Z (Q 2 ) (posAG). 



(2) 
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As(x, Q 2 )sidis differs essentially from the negative |(As + As)(x, Q 2 )bis ob- 
tained from all the QCD analyses of inclusive DIS data. It was also shown 
that when in the combined QCD analysis only the pion SIDIS data are in- 
cluded the polarized strange quark density is still negative for any x in the 
measured region, and the change-signing behavior of As (a;, Q 2 ) is due to the 
kaon asymmetries calculated by the DSS fragmentation functions. A further 
detailed analysis of the sensitivity of As(x, Q 2 ) to the kaon FFs is needed, 
and any model independent constraints on FFs would help. Another possible 
reason for this disagreement could be the assumption As(x, Ql) = As(x, Ql) 
made in the global analyses. However, this would probably require a sig- 
nificant difference between As and As, which is not seen in the COMPASS 
analysis. In any case, obtaining a final and unequivocal result for As(x) 
remains a challenge for further research on the internal spin structure of the 
nucleon. 

We have found also that the polarized gluon density is still ambiguous, 
and the present polarized DIS and SIDIS data cannot distinguish between 
the positive and a sign-changing gluon densities AG(x). This ambiguity is 
the main reason that the quark-gluon contribution into the total spin of the 
proton is still not well determined. 

Finally, our combined NLO QCD analysis confirms our previous results 
on the higher twist corrections to the nucleon spin structure function g± , 
namely, that they are not negligible in the pre-asymptotic region and have 
to be accounted for in order to extract correctly the polarized PDFs. 
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